Transport through single-wall metallic carbon nanotubes in the cotunneling regime 
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Using the real-time diagrammatic technique and taking into account both the sequential and 
cotunneling processes, we analyze the transport properties of single- wall metallic carbon nanotubes 
coupled to nonmagnetic and ferromagnetic leads in the full range of parameters. In particular, 
considering the two different shell filling schemes of the nanotubes, we discuss the behavior of 
the differential conductance, tunnel magnetoresistance and the shot noise. We show that in the 
Coulomb diamonds corresponding to even occupations, the shot noise becomes super-Poissonian 
due to bunching of fast tunneling processes resulting from the dynamical channel blockade, whereas 
in the other diamonds the noise is roughly Poissonian, in agreement with recent experiments. The 
tunnel magnetoresistance is very sensitive to the number of electrons in the nanotube and exhibits 
a distinctively different behavior depending on the shell filling sequence of the nanotube. 
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I. INTRODUCTION 

Transport properties of carbon nanotubes (CNTs) 
have been a subject of extensive studies for a few years ^ 
In particular, very recently there has been a growing 
interest in the shot noise measurements of CNTsj^ 5 ^ 
The shot noise provides useful information, not necessar- 
ily contained in the current, about the electronic struc- 
ture, couplings to external leads, and indicates the role of 
correlations and different types of processes driving the 
current. 7 In addition, recent experiments on CNTs cou- 
pled to ferromagnetic leads have also shown that nan- 
otubes exhibit a considerable tunnel magnetoresistance 
(TMR) effect^i^ii^^i^i^^l The TMR 
provides information about the spin accumulation on the 
nanotube and charge states taking part in transport ] 21 ! 22 

If the nanotube is weakly coupled to external leads, 
the current flows through the system due to consec- 
utive tunneling processes. The first-order (sequential) 
tunneling processes dominate above a certain thresh- 
old voltage and are suppressed below this voltage 
due to the single-electron charging energy and/or fi- 
nite level separation. On the other hand, current in 
the blockade regions is mainly mediated by second- 
order processes (cotunneling), 23 Theoretical considera- 
tions of spin-dependent transport properties of single- 
wall metallic CNTs in the perturbative regime have 
been so far mainly restricted to the sequential tunnel- 
ing processes i 2Q i 21 However, because transport in the 
Coulomb blockade regime is dominated by cotunneling 
processes, the sequential tunneling approximation may 
lead to wrong results. Therefore, to get reliable informa- 
tion on the transport properties in the full range of bias 
and gate voltages, which could be compared with that 
observed experimentally, one should go beyond the first- 
order theory. This is the main objective of the present 
paper. 

The considerations are based on the real-time diagram- 
matic technique, which allows us to take into account 



the sequential tunneling, cotunneling and cotunneling- 
assisted sequential tunneling processes in a fully sys- 
tematic way. Assuming realistic parameters of the 
system ^ 24 ! 25 we calculate the current, differential conduc- 
tance, TMR, and the shot noise of CNTs. In addition, we 
also analyze the effect of different shell filling sequences, 
which can be realized in the single-wall CNTs exhibit- 
ing four-electron periodicity ^ 24 ! 25 on transport character- 
istics. Generally, one can distinguish two different shell 
filling schemes - in the first scheme the following se- 
quence of the ground states is realized: S = 0, 1/2, 0, 1/2, 
where S is the spin of the nanotube, while in the second 
scheme it is: S = 0,1/2,1,1/2. We show that the ef- 
fect of different shell filling scenarios is the most visible 
in the TMR, while the current and shot noise are only 
slightly affected. This is due to the fact that the TMR 
directly reflects the magnetic properties of the system 
and it is thus very sensitive to whether the ground state 
of the nanotube is a singlet or a triplet. As concerns 
the shot noise, we show that the corresponding Fano 
factor F in the Coulomb blockade regions is enhanced 
above the Schottky value, F > 1, due to bunching of 
inelastic cotunneling processes. However, this enhance- 
ment is significantly lower than that obtained within the 
first-order approximation. 21 Contrary, in the transport 
regions where the sequential contribution to the current 
is dominant, the shot noise is sub-Poissonian with F 
slightly above 1/2 and is only weakly modified by the 
cotunneling processes. Our results for the shot noise are 
in qualitative agreement with recent experimental data 
on carbon nanotubes^ and quantum dots*2& 

The paper is organized as follows. In section 2 we 
describe the model of the nanotube and in section 3 
we briefly present the real-time diagrammatic technique 
used in calculations. Numerical results and their discus- 
sion is given in section 4, where we first present the re- 
sults on CNTs coupled to nonmagnetic leads and then we 
discuss the effects of spin-dependent tunneling on trans- 
port properties. In addition, we also present results for 
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two different shell filling scenarios. Finally, we give the 
conclusions in section 5. 



II. MODEL 

We consider a single wall metallic CNT weakly cou- 
pled to the two electrodes which can be either nonmag- 
netic or ferromagnetic. In the latter case the magne- 
tizations of the leads are assumed to be collinear and 
form either parallel or antiparallel magnetic configura- 
tion. The Hamiltonian H of the system takes the gen- 
eral form H = Hl + + Hqnt + The first two 
terms describe noninteracting itinerant electrons in the 
leads, H r = J^ka £rk<74k<7 C rk<7 for the left (r = L) and 
right (r = R) lead, with e r ^ a being the energy of an elec- 
tron with the wave vector k and spin a in the lead r. 
The third term of H describes the single-wall CNT and 
is given by2£ 

Hcnt = J2 £ M n M° + \ ( N ~ N of 

where N = E^v^wr, with n Mja = d^ ja d^ ja being the 
occupation operator for spin a and jth level in the sub- 
band /i (/i = 1,2). The energy e^j of the jth discrete 
level in the subband /i is given by e^j =jA+(/i — 
where A is the mean level spacing and S is the en- 
ergy mismatch between the two subbands of the nan- 
otube. The charging energy of CNT is described by U, 
and No is the charge on the nanotube induced by gate 
voltages. The additional on-level Coulomb energy be- 
tween two electrons occupying the same orbital level is 
denoted by 8U , whereas J describes the exchange en- 
ergy between the spin- up and spin-down electrons. The 
exchange effects described by J play an important role 
for small diameter nanotubes 2 ^ (up to 2 nm or so), 
as considered in this paper, whereas for nanotubes of 
larger diameters these effects become negligible. Fi- 
nally, the last term of H, Ht, takes into account tun- 
neling processes between the nanotube and electrodes, 

where t r j denotes the tunnel matrix elements between the 
lead r and the jth level (assumed to be spin-independent 
also for ferromagnetic leads). Coupling of the jth level to 
external leads is described by = 27r\t r j | 2 p^, with p a r 
being the density of states in the lead r for spin a. The 
role of ferromagnetic leads is taken into account just via 
the spin-dependent density of states p a r . The coupling pa- 
rameters can also be expressed as T^ - ^ = r rj -(l±p r ), for 
the spin-majority (spin-minority) electron bands, with 
T r j = (T+ +r~-)/2 and p r being the spin polarization of 
the lead r, p r = (/?+ — Pr)/(Pr + Pr )• For nonmagnetic 
leads, r^- = T~ . In the following we assume T r j = F/2 
for all values of j and r. 



III. METHOD 

In order to calculate the transport through a single- 
wall metallic carbon nanotube in the sequential and 
cotunneling regimes, we employ the real-time dia- 
grammatic technique j 21 i 29 i 3Q It consists in a systematic 
expansion of the nanotube (reduced) density matrix 
and the operators of interest with respect to the cou- 
pling strength T. The time evolution of the reduced 
density matrix can be visualized as a sequence of 
irreducible self-energy blocks, W, on the Keldysh 
contour. The matrix elements W xx > of W describe 
transitions between the many-body states |x) and |x') 
of the CNT*2i Then, the Dyson equation constitutes 
the full propagation of the reduced density matrix, 
which leads to a kinetic equation for the nanotube 
occupation probabilities, (Wp st ) x = T5 XXol where 
p st is the vector containing probabilities and the 
matrix W is given by W with one arbitrary row xo 
replaced by (T, . . . , V) due to the normalization con- 
dition ^2 x p x = 1. The current flowing through the 
system can be found from / = (e/2fi)Tr{W I p st }, 
whereas the zero-frequency current noise, S = 
2 dt((I(t)I(0) + J(0) /(*)> - 2</> 2 ), is given by^ S = 
(e 2 /h)Ti { [W n + W 1 (PW 1 + p st <8> e^W 1 )] p st }. 
The matrix W 1 ^ 11 ) is the self-energy matrix with 
one internal vertex (two internal vertices), resulting 
from the expansion of the tunneling Hamiltonian, 
replaced by the current operator, while the matri- 
ces <9W and dW l are partial derivatives of W and 
W 1 with respect to the convergence factor of the 
Laplace transform. 30 The object P is calculated from: 
WP = l(p st <g> e T - 1 - <9Wp st <g> e T ), with 1 being the 
unit vector with row xo set to zero, and e T = (1, . . . , 1) j^S 
In order to calculate the transport properties order by 
order in tunneling processes, we expand the self-energy 
matrices, = W^W 1 ) + W^X 2 ) + the 

dot occupations, p st = p st (°) + p st ( 1 ) + . . . , and, 
P = p( -1 ) -f p(°) + . . . , respectively. The self-energies 
can be calculated using the corresponding diagrammatic 
rules^^i 



IV. NUMERICAL RESULTS AND DISCUSSION 

The shell filling structure of single-wall carbon nan- 
otubes exhibits the four-electron periodicity, which is a 
consequence of the two spin-degenerate subbands of the 
nanotube. Furthermore, depending on the intrinsic pa- 
rameters of the nanotube, its shell filling can generally 
have two different scenarios ] 24 ! 25 The first one is associ- 
ated with the following sequence of the ground states: 
S = 0,1/2,0,1/2, which is realized when SU + J < S. 
On the other hand, in the second scenario the spin state 
of the nanotube changes as 5 = 0, 1/2, 1, 1/2, which hap- 
pens for SU + J > S. The difference between these two 
cases appears in the transport regime where the nanotube 
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FIG. 1: (Color online) Differential conductance (a) and Fano 
factor (b) as a function of bias and gate voltages. The pa- 
rameters are: A = 8.4 meV, 17/ A = 0.26, J/ A = 0.12, 
617/ A = 0.04, 5/ A = 0.27, k B T/A = 0.035, p L = Pr = 0, 
x = 0.14, and F = 0.2 meV. 



is doubly occupied. In the case of first sequence, i.e. for 
517 + J < 5, one orbital level of one of the subbands is 
fully occupied and the nearest orbital level of the other 
subband is empty - this is a singlet state. Whereas for 
517 + J > (5, there is a single electron on each level of the 
two subbands which, due to the exchange interaction J, 
lead to the formation of a triplet state. In the following, 
we analyze the behavior of the tunnel magnetoresistance 
and the shot noise in the linear and nonlinear response 
regimes depending on two different shell filling schemes. 



A. Shell filling sequence: S = 0, |, 



i 
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When 517 + J < 5, the shell filling of the nanotube 
is realized in this way that the next orbital levels are 
being occupied only after the lower lying levels are full, 
which leads to the sequence of doublet and singlet ground 
states, S = 0, 1/2, 0, 1/2. To model the nanotube in this 
transport regime, we have taken the experimental pa- 
rameters derived by W. Liang et al 24 Furthermore, to 
clearly elucidate the effect of cotunneling on transport 
properties, we first present and discuss the behavior of 
the differential conductance and the Fano factor in the 
case of CNT coupled to nonmagnetic leads. In this case, 
for comparison we also show the results obtained when 
taking into account only the first-order tunneling. 

In Fig. [1] we show the density plots of the differential 
conductance (a) and the Fano factor (b) as a function 
of the bias and gate voltages. Since the noise in the 
small bias regime (\eV\ < k^T), is dominated by ther- 
mal Nquist- Johnson noise, we excluded this part from 
considerations and marked it by the black thick line 
around V = (the corresponding Fano factor diverges as 
V — > 0). The black diamonds in Fig. [TJa) correspond to 
blockade regions, where sequential transport is exponen- 
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FIG. 2: (Color online) Bias dependence of the current, differ- 
ential conductance and the Fano factor for V g = —0.2666 V 
[(a)- (c)] and V g = -0.1202 V [(d)- (f)], see the dashed lines in 
Fig. [1] The dashed lines show the first-order calculation. The 
parameters are the same as in Fig. [T] 



tially suppressed and current is dominated by the second 
order processes. Furthermore, the four-electron periodic- 
ity is clearly visible in the linear conductance, while the 
additional lines in the differential conductance for volt- 
ages larger then the threshold for sequential tunneling 
are due to excited states participating in transport. 

To see clearly the contribution from cotunneling pro- 
cesses we show in Fig. [2] the current (a), conductance 
(b), and Fano factor (c) for two different values of the 
gate voltage, corresponding to vertical cross-sections in 
Fig. [1] through small and large Coulomb diamonds, see 
the dashed lines in Fig. [lj Additionally, in Fig. [2] we 
also show the corresponding quantities calculated in the 
first order (dashed lines). It is evident that the cotunnel- 
ing processes modify the current mainly in the blockade 
regime and also close to resonance (blockade threshold), 
leading to renormalization of the nanotube levels. This 
makes the conductance peak shifted towards smaller volt- 
ages. As regards the shot noise, cotunneling processes 
have a significant impact on the Fano factor, particu- 
larly in the blockade regions corresponding to even occu- 
pations, where the shot noise is super-Poissonian. The 
corresponding Fano factor is then larger than unity, but 
significantly reduced in comparison to that calculated in 
the first order approximation. 

Consider first the limit of sequential (first order) tun- 
neling processes. The sequential tunneling starts when 
the bias approaches the threshold voltage. Owing to the 
Fermi distribution, sequential processes start before the 
zero- temperature threshold is reached. When the high- 
est populated level is doubly occupied, tunneling from 
the source to the drain electrode may leave the nanotube 
either in the initial or in an excited state. When the 
CNT is left in the initial state, the next tunneling process 
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FIG. 3: (Color online) Tunnel magnetoresistance (a) and Fano 
factor in the parallel (b) and antiparallel (c) configuration as 
a function of bias and gate voltages. The parameters are as 
in Fig. Q]with p L = Pr = 0.5. 



goes with the same probability as the first one. When, 
however, the CNT is left in an excited state, the thermal 
(dynamical) blockade is lifted and the subsequent tunnel- 
ing processes are faster. This leads to bunching of fast 
tunneling processes, which leads to the super-Poissonian 



noise 
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The cotunneling processes, both single-barrier and 
two-barrier ones, are not blocked. Generally, there are 
elastic and inelastic ones, which occur with different 
probabilities. However, the asymmetry between the fast 
and slow ones is reduced in comparison with the first or- 
der, which leads to a smaller Fano factor when the second 
order processes are included. Apart from this, the inelas- 
tic cotunneling processes also lift the dynamical blockade 
of the sequential processes. 

It is easy to note, that the above scenario does not hold 
when the highest level is singly occupied. In such a case 
first order tunneling processes are of comparable velocity 
and therefore the Fano factor is only slightly above unity. 
This applies also to cotunneling processes. Thus, depend- 
ing on the ground state of the nanotube which can be 
changed by the gate voltage, we either find a considerable 
enhancement of the Fano factor (super-Poissonian shot 
noise) or the Fano factor is approximately equal unity 
(Poissonian shot noise). This observation is in agree- 
ment with recent experiments carried out by E. Onac et 
al£ Finally, it is noteworthy that outside the blockade 
regions the shot noise becomes sub-Poissonian and the 
corrections due to cotunneling processes are much smaller 
than in the blockade regime, although still noticeable. 
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FIG. 4: (Color online) Bias dependence of the current, differ- 
ential conductance, TMR and the Fano factor in the parallel 
(solid line) and antiparallel (dashed line) configuration for 
V g = -0.2666 V [(a)-(d)] and V g = -0.1202 V [(e)-(h)], see 
the dashed lines in Fig. [3] The parameters are the same as in 
Fig. El 



Consider now the numerical results on CNTs contacted 
to ferromagnetic leads, see Fig. OH The tunnel magnetore- 
sistance, Fig. [3^a), exhibits a distinctively different be- 
havior depending on the number of electrons on the nan- 
otube. In the linear response regime, between the two 
consecutive four-electron sequences, the TMR is given 
by the Julliere's formula, TMR = 2p L p R /(l — PlPr)* 37 
whereas in the other transport regimes it is suppressed 
below the Julliere's value. This is a completely new 
feature as compared to the sequential tunneling results 
where the linear response TMR is constant and given 
by TMR seq = plPr/0- ~ PlPr), irrespective of the gate 
voltage^Si 

In Fig. [3] (b) and (c) we show the density plots of the 
Fano factors in both parallel and antiparallel magnetic 
configurations. The general features of the Fano factors 
are similar to those in the limit of nonmagnetic leads, 
except for the fact that now there is an additional con- 
tribution to the shot noise coming from the difference 
between the two spin channels for transport. Because 
the resultant coupling of the nanotube to ferromagnetic 
leads in the parallel configuration is larger than in the 
antiparallel one, the fluctuations of the current are en- 
hanced in the parallel configuration as compared to the 
antiparallel one. This leads to the corresponding differ- 
ence in the Fano factors, see Fig. [3] (b) and (c), which is 
especially visible in the Coulomb blockade regime. We 
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FIG. 5: (Color online) Tunnel magnetoresistance (a) and Fano 
factor in the parallel (b) configuration as a function of bias 
and gate voltages. The parameters are the same as in Fig. [3] 
with 5/ A = 0.1. 



also note that although the magnitudes of the Fano fac- 
tors in both magnetic configurations are different, their 
general behavior is qualitatively similar. 

The bias dependence of the current, conductance, 
TMR and the Fano factor is shown in Fig. [H for two 
different values of the gate voltage and for both the par- 
allel and antiparallel magnetic configurations. All the 
quantities are calculated taking into account both the 
sequential and cotunneling processes, and we do not dis- 
tinguish between the first and second order contributions. 
As in the case of nonmagnetic leads, the second order 
processes significantly modify the TMR and Fano factor 
in the blockade regions, and only slightly outside these 
regions. An interesting feature is a drop of the TMR at 
the onset for sequential tunneling, see Fig. [H(c) and (g), 
which is present when the ground state of the nanotube 
is a singlet. This is due to the interplay of the inelastic 
cotunneling processes which flip the spin of electrons oc- 
cupying the nanotube and the sequential tunneling which 
is generally incoherent. On the other hand, when the 
ground state of the nanotube is a doublet, the TMR does 
not drop with approaching the threshold, see Fig. 03(a), 
which is due to a nonequilibrium spin accumulation in- 
duced in the nanotube. The spin accumulation tends 
to increase the TMR and, thus, reduces the role of the 
above-mentioned processes responsible for TMR drop in 
the case of ground state S = 0. 



B. Shell filling sequence: S 
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The intrinsic parameters of the nanotube depend 
mainly on their size and coupling to external leads i 24 i 25 
In particular, the quality of the coupling between the 
leads and nanotube contributes to the energy mismatch 
S between the two CNT's subbands. In the case when 
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FIG. 6: (Color online) Bias dependence of the current, differ- 
ential conductance, TMR and the Fano factor in the parallel 
(solid line) and antiparallel (dashed line) configuration for 
V g = -0.2544 V [(a)-(d)] and V g = -0.108 V [(e)-(h)], cor- 
responding to the dashed lines in Fig. [5] The parameters are 
the same as in Fig. [5] 



SU + J > 5, the sequence of the ground states is changed 
as compared to the case of SU + J < S. Now, in the 
charge state with two electrons in the nanotube, each 
electron occupies one orbital level of the two subbands. 
This, together with a finite J, leads to the formation of 
a triplet state in the nanotube. Thus, the shell filling 
sequence becomes S = 0, 1/2, 1, 1/2. 

In Fig. [5] we show the TMR and the Fano factor in 
the parallel configuration as a function of the bias and 
gate voltages. The Fano factor in the antiparallel config- 
uration is generally smaller than that in the parallel one, 
but its behavior is qualitatively similar, therefor we only 
show F in the parallel configuration - it is displayed in 
Fig. 03(b). Furthermore, it can be seen that the general 
behavior of the Fano factor is similar to that in the case 
of SU + J < S, see Figs. [U(b) and 02(b). Now, the Fano 
factor becomes super-Poissonian when the ground state 
of the nanotube is either a singlet or a triplet. This fact 
is rather intuitive as the dependence of the shot noise on 
the bias and gate voltages is rather conditioned by the 
charge states taking part in transport than by the mag- 
netic state of the system. The information about the spin 
state of the nanotube is mainly contained in the TMR, 
which for SU + J > S becomes considerably modified, es- 
pecially in the transport regime where the ground state 
is a triplet, see Fig. 03(a). 

The bias voltage dependence of the current, differential 
conductance, Fano factor in the parallel and antiparallel 



6 



configurations as well as of the TMR is shown in Fig. [6l 
One can now clearly see that general features of trans- 
port characteristics calculated for V g = —0.2544 V, i.e. 
when the ground state of the nanotube is a singlet, are 
similar to those shown in Fig. [U[a)-(d). The only differ- 
ence is the enhancement of the Fano factor at resonance, 
see Fig. f), and the drop of the TMR at the onset of se- 
quential tunneling where now TMR becomes slightly neg- 
ative, see Fig.[6jg). On the other hand, when the ground 
state of the nanotube is a triplet (V g = —0.108 V), the 
transport characteristics become even more modified, see 
Figs. Sfe) -(g) andE3(e)-(g). The main difference is visible 
in the behavior of the TMR which for SU + J > S is in- 
creased above the Julliere's value. This is associated with 
the nonequilibrium spin accumulation in triplet states, 
which tends to enhance the difference between the two 
magnetic configurations of the system. Furthermore, it 
is also due to a finite exchange interaction J, because of 
which only particular components of triplet take part in 
transport. In this way, the exchange interaction plays 
a role similar to an external magnetic field. As already 
shown in the case of quantum dots, the TMR may be 
enhanced above the Julliere's value when only particular 
spin states of the dot participate in transport. 34 

V. CONCLUDING REMARKS 

In conclusion, we have calculated transport character- 
istics, including shot noise and tunnel magnetoresistance, 
of CNTs weakly coupled to magnetic and nonmagnetic 
leads. We have also analyzed the effects of two different 



shell filling scenarios, which can be realized in single- wall 
metallic CNTs, on transport characteristics. In particu- 
lar, we have shown that the second order contributions 
are crucial in the blockade regions. Depending on the 
ground state of the nanotube, the shot noise is either 
supper-Poissonian due to the bunching of fast tunneling 
processes, or approximately given by the Poisson value, 
in agreement with recent experiments, 3 On the other 
hand, in the regions dominated by sequential transport 
processes the shot noise is sub-Poissonian with the cor- 
responding Fano factor roughly equal to 1/2. Further- 
more, we have also shown that the linear TMR exhibits 
a strong dependence on the gate voltage, being equal to 
the Julliere's TMR in transport regime between the two 
consecutive four-electron fillings and suppressed below 
that value in the other transport regimes. In addition, 
when the ground state of the nanotube is a triplet, the 
TMR in the nonlinear response regime at the onset for se- 
quential tunneling becomes enhanced above the Julliere's 
value, which is due to nonequilibrium spin accumulation 
in triplet states and finite exchange interaction J. 
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